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Abstract 

We investigate the effect of directed short and long range connections in a simple model of small 
world network. Our model is such that we can determine many quantities of interest by an exact 
analytical method. We calculate the function V(T), defined as the number of sites affected up 
to time T when a naive spreading process starts in the network. As opposed to shortcuts, the 
presence of un-favorable bonds has a negative effect on this quantity. Hence the spreading process 
may not be able to affect all the network. We define and calculate a quantity named the average 
size of accessible world in our model. The interplay of shortcuts, and un-favorable bonds on the 
small world properties is studied. 
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I. INTRODUCTION 



Real life networks, whether made by nature, (e.g. neural, metabolic and ecological net- 
works) , or made by human (e.g. the World Wide Web, power grids, transport 
networks and social networks of relations between individuals or institutes) jf| Q, Isj], have 
special features which is a blend of those of regular networks on the one hand and completely 
random ones on the other hand. To study any process in these networks, (the spreading of 
an epidemic in human society, a virus in the internet, or an electrical power failure in a large 
city, to name only a few), an understanding of their topological and connectivity properties 
is essential (for a review see [l| and references therein). Recently obtained data from many 
real networks show that like random networks hJEI; they have low diameter, and like reg- 
ular networks, they have high clustering. Since the pioneering work of Watts and Strogatz 



12j, these networks have attracted a lot of attentions and have been studied from various 
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In contrast to most of the models studied so far, many real networks like the World Wide 
Web, neural, power grids, metabolic and ecological networks have directed one-way links 
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19j,|20[. These types of networks may have significant differences in both their static 
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and dynamic properties with the Watts-Strogatz (WS) model and its variations 
The presence of directed links affects strongly many of the properties of a network. For 
example, for the same pattern of shortcuts, the average shortest path in an directed network 
is longer than that in an undirected one, due to the presence of bonds with the wrong direc- 
tions (blocks) in many paths. So is the spreading time of any dynamic effect on the lattice. 
Consider the quantity V(T) defined as the average number of sites which are visited at least 
once when we start a naive spreading process at a site and continue it for T steps. Note 
that we mean an average over an ensemble of networks and initially infected sites and by the 
naive spreading process we mean that at each step of the spreading process all the neighbors 
of an infected site are equally infected. The quantity V(T) may be taken as a crude ap- 
proximation for the number of people who have been infected by a contiguous disease after 
T time steps has elapsed since the first person has been infected. Clearly this is a simpli- 
fication of the real phenomena, since in real world a disease may not affect an immunized 
person or may not transmit with certainty in a contact. However as a first approximation, 
V(T) may give a sensible measure of the effect in the whole network. Since in an directed 
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network, an effect only spreads to those neighbors into which there are correctly directed 
links, there will be pronounced differences in this important quantity between an directed 
and an undirected network. As a concrete example consider a ring with N sites, without any 
shortcuts, where to emphasize the absence of shortcuts, we denote V(T) by Vq(T). If all the 
links have the same direction, we have Vo(T) = T, and if all of them are bidirectional, we 
have Vq(T) = 2T. In both cases the whole lattice gets infected after a finite time. However 
if the links are randomly directed then Vq(T) may be much lower and furthermore, there is 
a finite probability that only a small fraction of the whole lattice gets infected. 
Adding shortcuts to this ring of course has a positive effect on the spreading. In a sense 
we have a chance to see the interplay of two different concepts of small worlds in these 
networks. The size of the world as a whole may be small due to the ease of communication 
with the remote points provided by long range connections, however the world accessible to 
an individual may be small due to the absence of properly directed links to connect it to 
the outside world. 

It is therefore natural to ask how the presence of directed links and (or) directed short- 
cuts affects quantitatively the small world properties of a network? How we can make a 
simple model of a small world network with such random directions? A WS-type model 
for these networks may be as shown in figure However due to their complexity, these 
networks should usually be studied by numerical or simulation methods and they seldom 
amend themselves to exact analytical treatment. 



A. The aim, structure and results of this paper 

As we will show in this paper, with slight simplification one can introduce simpler models 
which although retain most of the small world features, are still amenable to analytical 
treatment. This is what we are trying to do in this paper. In this paper we introduce one 



such model following our earlier work 



2J. The 



23] which was in turn inspired by the work of 
basic simplifying feature of these networks is that all the shortcuts are made via a central 
site, figure (J2J). For such a network many of the small world quantities, can be calculated 
exactly. In particular, once V(T) defined above is calculated, many other quantities like the 
average shortest path between two sites can be obtained. An exact calculation of V(T) is 
however difficult for the case where both the shortcuts and the links have random directions. 
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FIG. 1: A WS-type model of directed network. 




FIG. 2: A simple substitute for the network of figure |[T|). 

We therefore proceed in two steps to separate the effects of randomness in the two types of 
connections. First, in section 2, we remove the shortcuts and calculate exactly V(T) for a 
ring with random links, figure (jHJ). To emphasize the absence of shortcuts we denote this 
quantity by Vo(T). Note that Vq{T) depends only on the structure of the underlying ring 
and its short-range connections. Then in section 3, we consider only the effects of randomly 
directed shortcuts, that is we let directions of the links on the ring to be regular and fixed 
say clockwise, and calculate exactly V(T), where again for emphasis on the shortcuts we 
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FIG. 3: A regular ring with randomly directed links, without shorcuts. You can see also the 
accessible world of site 1. 

denote this quantity by S(T). 

We then argue, in section 4 that in the scaling limit where the number of sites goes to 
infinity with the number of shortcuts kept finite, most of the spreading takes place via the 
links and only from time to time it propagates to remote points via the shortcuts. In this 
limit it is plausible to suggest a form for V(T) which takes into account the effect of both 
the random links and the shortcuts in the form V(T) = S(V (T)). This may not be an exact 
relation but as we will see it will give a fairly good approximation of V(T), as shown by the 
agreement of our analytical results and the results of simulations. This then means that in 
more complicated networks, one can separate the effects of short and long range connections 
and superimpose their effect in a suitable way. We conclude the paper with a discussion. 

II. EXACT CALCULATION OF V (T) IN A RING WITH RANDOM BONDS 

Consider a regular ring of N sites whose bonds are directed randomly. Each link may be 
directed clockwise with probability r, counterclockwise with probability £, and bidirectional 
with probability 1 — r — t. 

Thus we have a problem similar to bond percolation in a small world network. Suppose that 
at time T = site number 1 is infected with a virus. We ask the following questions: 
After T seconds how many sites have been infected on the average? What is the average 
speed of propagation of this decease in the network? These questions have obvious answers 
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for rings with regularly directed or bi-directional bonds, namely the number of infected sites 
are respectively T and 2T, with corresponding speeds of propagation being 1 and 2. In the 
randomly directed network, the situation is different. For example if both neighbors of site 
1 are directed into this site, this site can not affect any other site of the network. Such a 
Site be lng effect^ 1Sol a t e d has an m or U Q of Z e r o slZ e ( M To proceed 

with exact calculation, consider the right hand side of site 1. The probability that exactly 
k < T extra sites to the right have been infected is P+(k) := (1 — £) k £, and the probability 
that exactly T extra sites have been infected is P+{T) := (1 — £) T . Therefore the average 
number of extra sites infected to the right of the original site is 

V + (T) = Y.kP + (k)=T(l-£) T +J2k(l-£) k £ 

k=l k=0 

= T(l - £) T + j[l - £ + (1 - £f{£ - 1 - £T)\. (1) 
Going to the large iV limit where, 

T V + (T) 
N^oc, £^0 n:=£N, t := -, v + (t) := -2^i (2) 

we find the simple result 

v+{t) = kl-e-'*) (3) 
/i 

The same type of reasoning gives the number of sites infected to the left Vq (t) and thus the 
total number of infected sites will be: 

Vo (t) = -(l-e-^) + \(l-e- xt ) (4) 
/i A 

where A := rN. What are the meaning of the scaled variables? The parameter \i is the 
total number of sparse blocked sites in the way of propagation to the right, with a similar 
meaning for A. vo(t) is the fraction of infected sites up to time t. In a bidirectional lattice, 
all the sites could be infected after the passage of T = y seconds, or at t = | and if t 
passes |, some of the sites become doubly visited. Therefore it is plausible for the sake of 
comparison to define a quantity in our ring, namely the average size of the accessible world 
as Vq cc := fo(|), which turns out to be: 

v a cc = -(l-e-t) + l(l- e -i) (5) 
/i A 

It is seen that the presence of only a small number of blocked bonds causes a significant 
drop in the average size of this accessible world. For example a value of A = fx — 4 leads to 



Vq cc ~ 0.4. The long range connections (shortcuts) make the world small with the ease of 
communication they provide, however blockades make the world small in this new sense. 
The speed of propagation is found from 

v (t) = e-^ + e- xt . (6) 

In the symmetric case where A = fx, equation simplifies to: 

v (t) = -(l-e-n (7) 
A* 

with 

v Q (t) = 2e->* (8) 

Note that at the early stages of spreading when fit << 1, and the effects of blocked bonds 
has not yet been experienced, the infection propagates with speed equal to 2 as in a regular 
network. The effect of blocking comes into play when t becomes comparable to -. 
As a few number of shortcuts may enhance the speed of propagation, a few number of blocked 
bonds may have the opposing effect. First the blocks reduce the speed of propagation as is 
clear from (jUJ) and second and more importantly they reduce the number of accessible sites, 
or the size of accessible world. It will thus be of interest to see how these two effects compete 
in a random network where there are both shortcuts and blocks. We will study this in the 
final section of this paper. To this end we first study the effect of directed shortcuts in an 
otherwise regular ring with no blocks. 

III. THE LONG RANGE CONNECTIONS 

In this section we are to consider only the effect of randomly directed shortcuts in the 
spreading process and obtain exactly the function S(T) for this network, figure (jlj). Note 
that this function has the same meaning as V(T), except that for emphasis on the role of 
shortcuts in it we have adopted a new name for it. We fix a regular clockwise ring. Between 
a site and the center there is a shortcut going into the center with probability p and out of 
the center with probability q. The site remains unconnected to the center with probability 
1 — p — q. The average number of connections into and out of the center are respectively 
Mi := Np and M Q := Nq. 

Consider sites 1 and j. We want to find the probability that the shortest path between 
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FIG. 4: Randomly directed shortcuts added to a ring with clockwise links 

these two sites be of length I, a probability which we denote by P(l,j; I). A typical shortest 
path of length I connecting these two nodes is shown in figure (jlj, where the first inward 
connection to the center occurs at site i and the last outward connection from the center 
occurs at site j + i — l. Such a path occurs with probability (1 — p) t_1 pg(l — q) l ~\ Summing 
over all such configurations gives us the probability for the shortest path between sites 1 
and j to be of length I. For I ^ j — 1, we have: 

p(l,j; l± j - 1) = £(1 - P r l pq{l - q) 1 - 1 = pq[ { ^—^- + 2—^], (9) 
fri q-p p-q 

and p(l,j; j — 1) is determined from normalization: 

P(l,j; j - 1) = 1 - £ P(hf, I) = -^0(1 - Q) 3 ' 1 ~ 9(1 - P) 3 ~ l ) (10) 
z=i P ^ 

Note that p(l,j;l ^ j — 1) dose not depend on j, a property which is true for standard 
small world networks |25j ]. 

Now consider a naive spreading process starting at site 1. The number of sites affected 
up to time T, denoted by S(T), builds up in two ways, via the links on the ring and 
via the shortcuts. The first way gives a contribution T + 1 and the second way gives 
a contribution (N — T — 1) J2f=i p(l> j] 13 where (N — T — 1) is the number of sites 
beyond direct reach at time T which has been multiplied by the probability of any of these 
sites being at a distance shorter than T to site 1 via a shortcut. Putting this together we find: 
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FIG. 5: The speed of propagation in a ring, for several values of randomly directed shortcuts. 



T 



S(T) = T + 1 + (N-T-1)J2 P{hj; Q (11) 



i=i 

= N + (N — T— l)[-?-(l - p) T+1 + -£-(1 - q) T+1 ] 
p — q q — p 

In the scaling limit where N — > oo, p, g — > 0, where Mj and M Q are kept fixed and 

:= ^7^) we nnc l : 

s (t) = 1 _ ]~\ (Mie-™** - M e- Aht ) (12) 
In the symmetric case where Mj = M D = M this equation simplifies to: 

s(t) = 1- (1 - t)(l + Mt)e- Mt (13) 

with the speed of propagation 

s(t) = e~ Mt (l + Mt + M 2 t - M 2 t 2 ) (14) 
Figure (jSJ) shows the speed of propagation as a function of time for several values of M. 

IV. THE SPREADING EFFECT IN A DIRECTED SMALL WORLD NETWORK 

We now come to the problem of composing both the blocks and the shortcuts in a model 
of small world network. That is we consider the ring of figure (J2J) where randomly directed 
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shortcuts are added to a ring with randomly directed links. We can not obtain exact 
expressions for this network from first principle probability considerations. However we 
can obtain expressions for v(t) in the scaling limit by a heuristic argument and compare 
our results with those of simulations. Consider equation (|13|) . This equation shows how 
the presence of 2M randomly directed shortcuts in a regular clockwise ring affects the 
spreading effect. On the other hand we know that the number of sites infected up to time t 
in the absence of shortcuts, has changed to vo(t). Due to the rarity of shortcuts compared 
to the regular bonds, most of the spreading takes place via the local bonds, the role of 
shortcuts is just to make multiple spreading processes happen in different regions of the 
network. This role is the same whatever the underlying lattice is, and therefore for a general 
network, at least in the scaling regime, we can assume that equation (|13|) can be elevated 
to v(t) = s(v (t)), i.e; 

v (t) = 1 - (1 - v (t))(l + Mv (t))e- Mv °V . (15) 

For a fully random network where 2M randomly directed shortcuts are distributed on a ring 
with already random links, we assume that this relation holds true with Vo(t) taken from 
(0J). This suggestion may not provide an exact solution for the network, however we think it 
provides a fairly good approximation. In fact exact solution for the case where all the links 
on the ring are bidirectional is possible and it confirms the above ansatz, that is we obtain an 
exact expression only by setting Vo(t) = 2t in the above formula. Moreover this separation 
of the effect of short and long range connections may be also useful in more complicated 
networks. Whether this assumption is plausible or not can be checked by comparison with 
simulations. The results of simulations are compared with those of equations (J3]) and (|T2|) in 
figure USD and (0). 

V. STATIC PROPERTIES 

Once the functions V(T) or v(t) are obtained, the static properties of the network i.e., 
the average shortest path between two arbitrary sites and its probability distribution can 
be calculated directly. 

Since V(T) by definition is the number of sites whose shortest distance to site 1 is less 
than or equal to T, we find the number of sites whose shortest distance is exactly T to 
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FIG. 6: V(T) for a fully random network in the case N = 5000, r = 0.02,1 = 0. Analytic 
results (lines) versus simulations(symbols) which have been averaged over 1000 realizations of the 
network. 
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FIG. 7: V(T) for a fully random network in the case N = 5000, r = I = 0.005. Analytic re- 
sults(lines) versus simulations(symbols) . 

be V(T) — V(T — 1). Since site 1 is an arbitrary site, we find the probability distribution 
of the shortest distance between two arbitrary sites which are accessible to each other as: 
P(T) = YJUzYVjz^) w here V acc is the average size of the accessible world. (There is of 
course a slight approximation here in that we are taking averages of the denominator and 
numerator separately.) 

For a regular ring with shortcuts, V acc = N, since all the sites are accessible. We will discuss 
the case of random rings in the sequel. In the scaling regime the above formulas transform 
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to: 

P{t) = v(t). (16) 

Note that P(t) is normalized, i.e. /„ P(t)dt = v(l) - w(0) = 1. The average shortest path 
for the network of figure (j3J) when M, = M = M, turns out to be: 

(t) = f 1 t P(t)dt = f 1 tv(t) = ^-(2M-3 + (M + 3)e- M ) (17) 



This is in accord with the result of |24J]. This formula shows that the presence of a small 
number of shortcuts, causes a significant drop in the average shortest path from 1 to very 
small values. In this sense the world gets smaller by long range connections. 
We now study the static effects of random directed bonds on a ring without shortcuts. The 
presence of blocks makes the world small in a different sense, namely for each site the number 
of accessible sites gets smaller. In fact the average size of the world accessible to a site is 
not N anymore but it is given by V(y) (see the paragraph leading to equation ©). Hence 
the probability of shortest paths is given by P{T) := v ^~Yn\~ 1 \ or in the scaling limit by 

Pit) := 3 (18) 

1 

This probability is normalized, i.e. J 2 P(t)dt = 1. We obtain from (j!8)l 



(0 = ~~ 7TT / tv(t)dt (19) 
v{±) Jo 

However in order to assess the situation in this network, we should compare the average 

shortest path with the size of this small world itself, namely we should calculate In- 
serting equation (|ZJ) into (fTHjl we find: 

= 2 -<" + 2 ^ (20) 

vS cc 4(1 -e-f) 2 V ; 

Figure (jHJ shows both the average size of the accessible world Vq cc and the ratio of the 
average shortest path to the size of accessible world as a function of the number of blocks \i. 
It is seen that for /j, — 0, when there is no block, the size is 1 and the average of the shortest 
path is I as it should be. With a few number of blocks the size drops dramatically and the 
average of shortest path within the world increases. Note that with increasing jj, the average 
shortest path increases to its maximum value of ^. 

For the fully random network, we use equations (JT5j) and (|18|) to obtain the average of 
shortest path. The result is shown in figure © for several values of the parameters. 
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FIG. 8: The average size of accessible world and the average shortest path for a regular ring with 
randomly directed bonds without shortcuts. 
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FIG. 9: The average shortest path for a fully random network. 
VI. CONCLUSION 



We have studied the effect of directed short and long range connections in a simple model 
of small world network. In our models all the shortcuts pass via a central site in the network. 
This makes possible an almost exact calculation of many of the properties of the network. 
We have calculated the function V(T), defined as the number of sites affected up to time 
T when a naive spreading process starts in the network. As opposed to shortcuts, the 
presence of un-favorable bonds has a negative effect on this quantity. Hence the spreading 
process may be able to affect only a fraction of the total sites of the network. We have 
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defined this fraction to be the average size of the accessible world in our model and have 
calculated it exactly for our model. We have studied also the interplay of shortcuts, and 
un-favorable bonds on the small world properties like the size of accessible world, the speed 
of propagation of a spreading process, and the average shortest path between two arbitrary 
sites. Our results show that one can separately take into account the effect of randomness 
in the directions of shortcuts and the short-range connections in the underlying lattice and 
at the end super-impose the two effects in a suitable way. We expect that this will hold also 
in more complicated lattices of small world networks. 
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